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A black hole has characteristic quasi-normal modes that will be excited when it is formed or when the geometry 
is perturbed. The state of a black hole when the quasi-normal modes are excited is called the gravitational ringing, 
J> ' and detections of it will be a direct confirmation of the existence of black holes. To detect it, a method based on 
\ matched filtering needs to be developed. Generically, matched filtering requires a large number of templates, because 
one has to ensure a proper match of a real gravitational wave with one of template waveforms to keep the detection 
efficiency as high as possible. On the other hand, the number of templates must be kept as small as possible under 
. limited computational costs. In our previous paper, assuming that the gravitational ringing is dominated by the 
' least-damped (fundamental) mode with the least imaginary part of frequency, we constructed an efficient method for 
tiling the template space. However, the dependence of the template space metric on the initial phase of a wave was 
' not taken into account. This dependence arises because of an unavoidable mismatch between the parameters of a 



^H, signal waveform and those given discretely in the template space. In this paper, we properly take this dependence 



^ ■ into account and present an improved, efficient search method for gravitational ringing of black holes. 



I. INTRODUCTION 



5-H ' Thanks to the recent technological advance, there are many on-going projects of gravitational wave detection in the 



world; the Laser Interferometric Gravitational Wave Observatory (LIGO) [1], VIRGO [2], GEO-600 [3], ACIGA [4], 
TAMA300 [5] and the Large-scale Cryogenic Gravitational wave Telescope (LCGT) [6] which are ground-based laser 
interferometers, and EXPLORER [7], ALLEGRO [8], NIOBE [9], NAUTILUS [10] and AURIGA [11] which are 
bar detectors. Furthermore, there are some future space interferometer projects such as the Laser Interferometer 
Space Antenna (LISA) [12] and the DECi hertz Interferometer Gravitational wave Observatory (DECIGO) [13]. The 
detection of gravitational waves provides us with not only a direct experimental test of general relativity but also a 
new window to observe our Universe. To use them as a new tool of observation, it is necessary to derive theoretical 
waveforms. Once we know them, we may appeal to the matched filtering technique to extract source's information 
from gravitational wave signals. However, because the signals are expected to be very weak and the amount of data 
will be enormous for long-term continuous observations, it is essentially important to develop efficient data analysis 
methods. 

For these ground-based as well as future space-based interferometers, coalescence of compact object binaries is 
the most important source of gravitational waves. The process of a binary coalescence can be divided into three 
distinct phases. During the initial inspiral phase, the gravitational radiation reaction timescale is much longer than 
the orbital period. The gravitational waves from the inspiral phase carry the information of the masses, spins and so 
on, of the system. After the inspiral phase, the binary becomes dynamically unstable and starts to merge. This phase 
is called the merger phase. The gravitational waves from the merger phase give us the information about fully general 
relativistic dynamics of the system. Finally, if a black hole is formed after the merger phase, the system enters the 
ringdown phase where it gradually settles down to a stationary Kerr black hole. During this process, the black hole 
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emits gravitational waves with frequencies and damping rates specific to its mass and spin. Thus, the gravitational 
waves in the ringdown phase carry the information of the mass and spin of the final black hole. 

In this paper, we consider an effective search method for gravitational ringing of distorted spinning (Kerr) black 
holes. The ringdown waves are described by quasi-normal modes of a black hole. The quasi-normal modes are complex 
frequency wave solutions of the perturbed Einstein equations with purely outgoing- wave boundary condition at infinity 
and ingoing-wave at horizon, with vanishing incoming-wave amplitude. A quasi-normal mode is characterized by the 
central frequency fc, usually called the (quasi-)normal-mode frequency, and the quality factor Q which is inversely 
proportional to the imaginary part of the complex frequency. Their properties were analyzed extensively by Leaver [14], 
and it is known that the least-damped (fundamental) mode belongs to the £ = m ^ 2 spin-2 spheroidal harmonic 
modes, we note that a ringdown signal decays exponentially. Therefore, unless the signal-to-noise ratio (SNR) is very 
large, the signal will be soon buried in the noise after a few oscillation cycles. So, it is essentially important to keep 
the loss of SNR as small as possible when wc construct a search method. 

Here we consider a black hole characterized by its mass M and the dimensionless spin parameter a = J /M^ where 
J is the spin angular momentum. The parameter a takes a value in the range [0, 1), with a = corresponding to 
a Schwarzschild black hole and a = 1 to an extreme Kerr black hole. It is known that the imaginary part of the 
i = m = 2 least-damped mode is the smallest of all the quasi-normal modes, and results of black hole perturbation 
calculations as well as numerical relativity simulations strongly suggest that a ringdown wave is dominated by this 
£ = m = 2 least-damped mode unless a is extremely close to unity. Hence we focus on this single mode. Then a 
ringdown waveform is expressed as 

Kfc, Q, to, (t>o; t) = [ e" ^"'^ cos(2 7r/c {t - to) - 0o) for t > to , 
I for t < io , 

where we set the amplitude to unity for simplicity, and to and (po are the initial time and phase of the ringdown 
wave, respectively. For the £ = m = 2 least-damped mode, analytical fitting formulas for the central frequency fc and 
quality factor Q were given by Echeverria [15] as 

/e ~ 32kHz [1 - 0.63 (1 - af^] (^^) , (1.2) 
Q ~ 2.0(l-a)-°•*^ (1.3) 

There exists rich literature on search methods for ringdown waves. Echeverria [15] investigated the problem of 
extracting the black hole parameters from gravitational wave data in the case when SNR is large. Finn [16] improved 
this situation by developing a maximum likelihood analysis method that can deal with any SNR. Flanagan and Hughes 
then considered the parameter extraction from the three stages of a binary coalescence, i.e., from inspiral, merger 
and ringdown phases, in their series of papers [17]. For the ringdown phase, they discussed the relation between 
the energy spectrum of the radiation and SNR. Creighton [18] reported the result of analyzing data of the Caltech 
40m by matched filtering, and emphasized the importance of coincidence event searches to discriminate spurious 
events from real events. But the search was limited to a single ringdown wave template. Recently, Arnaud et al. [19] 
discussed a tiling method to cover the 2-dimensional template space {fc, Q}- In our previous paper [20], we proposed 
a more efficient method for tiling the template space. There, however, we ignored the dependence of the metric of 
the template space {fc, Q} on the initial phase (f)Q. This induces some small but non-negligible decrease in the match 
for a signal with certain ranges of ^o- In this paper, we remove this shortcoming by properly taking into account 
the initial phase dependence, and develop a similar but substantially improved template spacing which is much more 
reliable than the previous one. 

Here we make a comment on an analysis using real interferometers' data. In this case, we have to deal with non- 
stationary, non-Gaussian noises, and Tsunesada and Kanda [21] found that more fake events are observed than in the 
case of an inspiraling wave search, because the duration of a ringdown wave is typically much shorter than that of 
an inspiral wave and it can easily be affected by short bursts. It will be necessary develop a way to remove such fake 
events without losing real ringdown signals, but we leave this issue for future work. 

The paper is organized as follows. In Sec. II, first, we introduce orthonormal template waveforms for ringdown 
waves. Second, assuming a white noise background, we consider matched filtering in the 4-dimensional template space 
{^0, ^Oj fc, Q}- We show that we can effectively reduce the template space to 2-dimensions spanned by {fc, Q}, but 
the metric of this reduced template space depends on <po- Then, by carefully taking account of the (f)o dependence 
of the metric, wc analytically develop an efficient and reliable tiling method for ringdown wave searches. In Sec. Ill, 
by using a fitting curve for the TAMA noise spectrum during the Data Taking 8 (DT8) in 2003, we show that our 
template spacing developed for white noise is valid even in the case of colored noise. Finally, Sec. IV is devoted to 
summary and discussion. In Appendix A, we recapitulate the tiling method we proposed in [20]. In Appendix B, we 
summarize the parameter estimation errors for ringdown signals by using the Fisher information matrix. 
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II. TEMPLATE SPACE 

In this section, we present an improved version of the template spacing we have developed in our previous paper [20] , 
which can be used for matched filtering of the quasi-normal ringing waveforms. Here, the detector noise is assumed 
to be white noise to make it possible to deal with the problem analytically. 

A. Maximization over the initial phase 00 

We have temporarily set the amplitude to imity for simplicity in Eq. (1.1). Note that the knowledge of the amplitude 
is not necessary for the template spacing in matched filtering. For t >to, the ringdown wave (1.1) is divided into two 
parts. 

h{fc, Q, to, 4>o; t) = hcifc, Q, to; t) cos^o + hs{fc, Q, to; t) sm4>o , (2.1) 

where 

/ic(/c, Q, io; t) = e" cos(2 7r/e (i - io)) , (2.2) 
hsifc, Q, to; t) = e-"^°«"*°' sin(2 Trf^t- to)) . (2.3) 

Performing the Fourier transformation haif) = dte'^'^^^*h„{t), where cr = c or s, we obtain the waveforms in 
the frequency domain as 

h(f Off)- (/c-2»/Q)ge^-/«" 

ncUc,Q,to,n - ^(2/^Q_,/^_2/0)(2/,Q + z/e + 2/Q)' ^''^^ 

hsifc, Q, to; f) = ^(2/,Q-i/,-2/g)(2/eQ + i/e + 2/g) • ^^-^^ 
The waveform in the time domain is real, so the following relation is satisfied. 

h*{f)=h{-f), (2.6) 

where the star ( *) denotes the complex conjugation. 
Here, we introduce the inner product as 

where Sn is defined as a one-sided power spectrum density, and /max is the maximum frequency we take into account 
in the analysis. In the actual data analysis, it is equal to or less than the half of the sampling frequency of data. In 
this section, the detector noise is assumed to be white noise Sn{\f\) = 1- We will use the fitting curve for the TAMA 
DT8 noise spectrum as a colored noise in the next section. 

In the matched filtering, we calculate the inner product between the template h and a signal x defined by {x,h). 
First, we normalize the templates he and hg. We define the normalization constants as 



Ncifc, Q, to) = {hcifc, Q, to), hcifc, Q, to)) 
^ (2Q^ + 1)Q 
7r(4Q2 + i)/^' 

Nsifc,Q,to) = ihsifc,Q,to),Kifc,Q,to)) 
^ 2Q3 

7r(4Q2 + i)/^' 



(2.8) 

(2.9) 



when we have set /max = oo. It is noted that the normalization constant does not depend on the initial time. The 
normalized templates h^ifc, Q, to) are given by 

Life, Q, to; f) = -^==^==Kifc, Q, to; f) • (2.10) 

V J^aifc, Q, to) 
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We note that the two parts, he and hg, are not orthogonal. Their inner product is obtained as 

(hcifc, Q, to),hsifc, Q, to)) = 



\/2(2Q2 + i) 
=: c(/„ Q, io) . (2.11) 

Then the normahzed ringdown template for the waveform (1.1) is given in the frequency domain as 

k^o; fc, Q, to; f) = -y==l===h{fc, Q, to, 0o; /) , (2.12) 

\/N{(j)o; Jc, Q, to) 

where 

h(f n f ^ f\ - (cos (/>o /c - 2 z cos 00 / Q + 2 Q sin 0o) g ^ " *° 
ftUc, W, to, 00, /j - ^(2/,Q-i/,-2/Q)(2/,Q + i/, + 2/Q) ' 

-^''(';^o; fc, Q, to) = {h{fc, Q, to, 0o),/i(/c, Q, to, M) 

= Ncifc, Q, to) cos^ 00 + Nsifc, Q, to) sin^ 0o 

+2c(/c, Q, to)VNc{fc, Q, to)Ns{fc, Q, to) cos 00 sin 00 . (2.13) 

In an actual data analysis, instead of the normalized templates he and hg, it is more convenient to prepare a set of 
orthonormalized waveforms hi and /12 as templates. They can be obtained by Schmidt's orthonormalization. First, 
we choose 

hi{fo,Q,to) = k{fc,Q,to). (2.14) 
Then the normalized waveform /i2 orthogonal to hi is obtained as 

, If hs{fc,Q,to)-c{fc,Q,to)hc{fc,Q,to) 

VI - c(/c, Q, to)^ 

Using the above orthonormal templates, the SNR maximized over the initial phase 0o is given as 

p{x;fc, Q, to) = Max|^(a;,/i(/c, Q, to, 0o))2| = {x, hi{fc, Q, to))^ + {x, h^ih, Q, ^o))^ (2.16) 

where x represents a signal. We note that there is basically nothing wrong with using the non-orthogonalizcd templates 
ha- in the matched filtering. In the case we use the non-orthogonal templates, the maximization of {x, h) over the 
phase 00 has been discussed by Mohanty [22], and the result is the same if we replace hi and ft,2 in Eq. (2.16) by 
those given in Eqs. (2.14) and (2.15), respectively. 

Before closing this subsection, let us describe how to determine the initial phase 0o from the filtered output. For 
simplicity, we take the normalized wave /i as a signal, because the overall amplitude is irrelevant for the present 
discussion. We further assume that there is a perfect match of the parameters {fc, Q, to} between the signal and a 
template, and the SNR is maximized, that is, p{h; fc, Q, to) = 1- 

We consider the quantities, 

pi = {h, hi) , 

P2 = {k h2) . (2.17) 
When the signal perfectly matches with a template, these are formally expressed as 

pi = ^ = (VNcifc, Q, to)cos0o + c(/e, Q, to)V^.(/c, Q, to)sin0o) (2.18) 

V^i'po; fc, Q, to) ^ ' 

1 



P2 



V7V(0o; /c, 0, to) 
Taking the ratio of pi and p2, we find 



v/(l - c(/c, Q, to)^)Ns{fc, Q, to) sin 00 . (2.19) 



Bl= / Nc{fc,Q,to) / c{fc,Q,to) 

P2 V (1 - c(/c, Q, to)2) N,{fc, Q, to) ^° ^ V 1 - c(/c, Q, to)2 ■ ^ ' ' 



5 



Thus, the initial phase of a ringdown wave signal is determined as 

4>o = cot~^ 



I{1 - c{fc, Q, toy)Ns{fc, Q, to) pi c{fe, Q, to)Ns{fc, Q, h) 



Ncifc, Q, to) 
It is noted that (^o ^ cot~^{pi/p2) for large Q. 



P2 



Ncifc, Q, to) 



(2.21) 



B. 3-dimensional distance function 

In the previous subsection, we obtained the expression for the SNR maximized over the initial phase (po, Eq- (2.16). 
Here, we consider the remaining 3 parameters {/c, Q, to}, and derive the metric in the 3-dimensional template space 
that describes the degree of mismatch between a signal and a template. To this end, we point out an important fact 
that was overlooked in [20]. Given a signal with certain values of the parameters, the SNR defined by Eq. (2.16) will 
be independent of (f)o for a template that matches exactly with the signal. However, it will not be so if there is no 
exact matching between the signal and a template. In other words, if there is a mismatch between the parameters of 
a signal and those of a template, the resulting SNR will depend on the initial phase 0o- 

Let us define the match C(0o; dfc dQ, dto) between the template with the parameters (/c, Q, to, 4>o) and the 
normalized signal having slightly different values of the parameters (/c + dfc, Q + dQ, to + dto, (l>o)- Note that, since 
we have already maximized the SNR over the initial phase (po, it is unnecessary to consider the difference dcpo for the 
match. The match C(^o; dfc, dQ, dto) is defined by 



C((^o; dfc, dQ, dto) 



P 



+ dfc, Q + dQ, to + dto, M\ fc, Q, to ' 
y-v \h{fc + dfc, Q + dQ, to + dto, (po), hA{fc, Q, to))^ 



[ {Kfc + dfc, Q + dQ, to + dto, <Po), Kfc + dfc, Q + dQ, to + dto, <Po)) 



1/2 



{h,i{fc, Q, to, (po), hj{fc, Q, to, (po)) 



2(/i(/c, Q, to, (po), Kfc, Q, to, (po)) . 

2 

- ^{h,i{fc, Q, to, (po), hAifc, Q, to)){hj{fc, Q, to, (po), hAifc, Q, to)) 

A=l 

+0{dx^) 

1 - ds^3) + 0{dx^) , 



dx^dx^ 



(2.22) 



where {x'} = {to, fc, Q} (* — 1,2,3), the comma (,) denotes the partial differentiation, and we have introduced the 
3-dimcnsional distance function in the template space by ds"^^^ = 1 — C. The smaller the match C is, the larger the 
distance is between the two signals in the template space. 
If we introduce the metric [23] as 



d4x) =g^fdx'dx\ 
the components of the metric are explicitly written as 

(4Q2 + 1) COs2(0o) 



(2.23) 



9toQ 



Q (4 Q2 + 2 g sin(2 </)o) + 1 + cos(2 0o)) 



TT (4Q2 + 1) sin(2 0o) 



2 (4 Q2 + 2 Q sin(2 (Po) + l + cos(2 .^o)) 
Ho 

IT fc cos(0o) (2Q cos((/>o) - sin(</>o)) 



(3) 
9Qto 



Q (4 Q2 + 2 Q sin(2 <^o) + 1 + cos(2 ,^o)) ' 
g2 (4 Q2 _ 2 Q sin(2 0o) + 1 - cos(2 4>o)) 
2 fc (4 Q2 + 2 Q sin(2 <^o) + 1 + cos(2 ,^o)) ' 
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„(3) _ _(3) 

9fcQ - 9Qf, 



Q{2Q cos(0o) - sin((/)o))^ 



(3) ^ K--^ ' ' V ^ — v-Tuy . - -t. v--t -/ y-T^; ,^ 24) 

2(4 02 + 1)2 (4 Q2 + 20 sin(2 0o) + l + cos(20o)) 

where we have assumed /max ^ /c and retained only the leading terms in /max//c- Note that /max appears only in 
the component g^^l^ in this limit. In the next subsection, we shall argue that the dependence on /max will disappear 
from the final result if the maximization over to is done with accuracy dto ^ l//max- 

The inequality C(0o; dfc, dQ, dto) < 1 means that there will be a loss of SNR unless the actual parameters of a 
gravitational wave signal fall exactly onto one of the templates. In order to keep the detection efficiency high enough, 
we have to choose the template spacing so that the maximum value of the distance function is kept as small as possible. 



fc (4 0" + 1) (4 0' + 2 sin(2 (j)o) + I + cos(2 (j)o j) ' 
{AQ^ + lf + (12 0^ - 1) cos(2 4>o) + 2Q (4 - 3) sin(2 ^o) 



C. Projection to 2-dimensional template space 

The 3-dimensional metric on the template space {to, fc, Q} has a simple /max dependence in the large /max limit, 
namely, gl^l^^ cx /max and the other components arc independent of /max- This fact allows us to reduce the problem 
of template spacing to that in the 2-dimensional space {fc, Q}- 

To reduce the 3-dimensional template space to 2-dimensions, let us note the distance function dsj,^^ may be rewritten 
as [23] 

= gill (dto + ^dx'] + gijdx'dx' , (2.25) 



9toto 



where {x^} = {fc, 0} and 



„(3) (3) 

9IJ = gf] - -^^^ ■ (2.26) 

We see that we can minimize the distance function (maximize the match) with respect to to by filtering a data stream 
with respect to to with accuracy dto = 0(l//max), because we then have 

gill [dto + ^dx'\ = 0(/c//max) ^ , (2.27) 

in the limit fma-K/ fc 3> 1. Thus, provided the matched filtering is done with this accuracy, which is normally practiced 
when we sweep over a data stream in time sequence, we may focus on the 2-dimensional space {fc, 0} for the problem 
of template spacing. Note that in the limit fma^/fc — * oo, the reduced 2-dimensional distance function is given by 

ds^2) = gijdx^dx-' = g'fjdx^dx'^ . (2.28) 

Thus the /max dependence disappears completely. 

Before proceeding further, we note that the metric gu has a simple central frequency dependence which may be 
removed by the coordinate transformation, 

dF = |; 
fc ^ F = ln(/e//o), 

F ^ fc = foe^ , (2.29) 
where fo is some fiducial frequency. This transformation gives 

rfs(2) = gPF dF^ + 2 gpQ dF dQ + gqq dQ^ ; 

gPF = flgu,, 9FQ=fc9f,Q- (2.30) 
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Now, to find an appropriate tiling of the 2-dimensional template space, we consider contours of a fixed maximum 
distance ds^^) = '^■^max- ^^ote that the metric has the initial phase dependence. So, the contour varies as (j)o is varied 

(over the range < < tt). The contours for ds'^^^ = 0.02 are plotted in Fig. 1 in the case of the quality factor 
Q = 2. The thick curve shows the contour with (j)o = which we considered in [20]. We sec that if wc use the previous 
template spacing given there, the minimum match is not guaranteed for other values of the initial phase (/>o- It is 
therefore necessary to remedy this drawback when tiling the template space. 




FIG. 1: Contours of the maximum distance dSmax = 0.02 for various initial phases. The quality factor is set to Q = 2. The 
thick curve shows the case of (po = 0, the interior of which corresponds to the region covered by one template in [20]. 



D. 2-dimensional template spacing 

We now consider the template spacing that guarantees the minimum match for any value of 4>o- As seen from 
Fig. 1, a conservative method is to take the interior of the inner envelope curve of the contours to be a region covered 
by each template. If wc adopt such a template spacing based on the exact inner envelope curve, the computational 
cost of assigning grid points in the template space becomes high. We therefore choose to allow some redundancy, and 
construct an ellipse that can cover the interior of the inner envelope curve as much as possible. This also enables us 
to use the efficient method of template space tiling developed in [20] . 

Let us determine the inner envelope curve. It is convenient to perform a coordinate transformation that makes the 
envelope curve symmetric with respect the axes of new coordinates. To do so, we first make a rotation in the {F, Q)- 
plane to bring the two points of intersection of all the contours on to the axes. This is achieved by the coordinate 
transformation 

du = cos 6 dF + sin 9 dQ , 

dv = - sine dF + COS0 dQ , (2.31) 

where 

1 

cos^ : 



The metric components in the new coordinates are given by 

^8^2) = 9uu du^ + 2 Quv du dv + g^v dv"^ ; 



sinO=^^M±^=. (2.32) 
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Quu = cos^ 9 gpF + 2 cos 6* sin gpQ + sin^ 9 gqq , 
9uv = — cos9 sm9 gpp + {cos^ 9 ~ sin^ 9) gpQ + cos6 sm9 
gyy = SVC? 9 gpp — 2 cos9 sm9 gpQ -\- cos^ 9 gqq . (2.33) 

Using the large Q expansion valid for Q >2, the new coordinates u and v as functions of Q, to 0{1/Q^) inclusive, 
are given by 



1 1 
+ 



4 legs 64 g7 



F 



Q + 



1 



1 



V = — 



32 06 64 0^. 
The inverse transformation is given by 

1 1 



160 0^ 448 0^ 
1 1 1 



+ 



8 0=^ 64 0* 384 0^ 2048 0^ 



F = 



Q = - 



-1 



32 «6 64 u8 



1 



1 



1 



+ 



39 



16 u7 



1 



4u3 16 m5 64 u7 



1 1 
+ 



8u2 64 m4 384^6 10240^8 
1 17 



u + 



40 u5 1792^7 



The contours of the maximum distance ds-^^ = 0.02 for = 2 on the {u, i')-plane are shown in Fig. 2. 



(2.34) 



(2.35) 



-.iM 



-O.S -0.4 -0.2 



FIG. 2: Contours of the fixed maximum distance dSmax = 0.02 for Q = 2 on the (w, v)-plane with the metric (2.33). 



Next, we derive the envelope curves. An envelope curve is the solution of the following two equations. 

R{du, dv; 0o) = > 



d 



R{du, dv; 4>q) = , 



(2.36) 



where the function R is given in our situation by 

R{du, dv; (po) — guu dv? + 2 g^v du dv + g^y dv^ — '^■^max ■ 
The envelope curve is obtained by eliminating the initial phase (po from Eqs. (2.36). We find 



(2.37) 



du + 



Q(4Q2 + 1) 



dv 



16 0*^ + 80^ + 0^ + 1 
4 0^(4 02 + 1)2 
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2Q(16g^ + 8(52 + l)rfsLx-<3(16Q^ + 8g^ + g2 + l)rfv^ 
±(4 g2 + 1) V2 (16 Q6 + 8 Q4 + Q2 + 1) dg^^ i^^i ^ 



(2.38) 



where Q is a function of u and u as given in Eq. (2.35). The outer and inner envelope curves are given by the plus and 
minus signs, respectively, in front of the square root in the above solution. Since our interest is in the inner envelope 
curve, we focus on the minus sign solution. 

As we can see from Eq. (2.38), the envelope curve is still asymmetric with respect to both axes. To make it 
symmetric, we make a further coordinate transformation. 



dV = dv, 
dU = du-\ 



'(4g2 + i) 



dv . 



(2.39) 



which gives the inner envelope curve in the form, 

16Q6 + 8Q4 + q2_^^ 



2Q{16Q^ + 8Q^ + l)ds 



2 

max 



QiUQ"" + 8Q^ + + l)dV^ 



-(4Q2 + 1) ^/2^T6Q^T8Q^TQ^TT)ds„,^^\dV\ 



(2.40) 



The equations (2.39) are integrated to give 
V = V 



1 - 



U 



Q 



32 

1 



+ 



1 



+ 



64 Q8 
1 



F- 



1 1 
+ 



1 



+ 



8Q2 64 Q4 384 Q6 2048 



160 448 
The inverse transformation is given by 

1 1 



(2.41) 



F = 



-1- 



Q = U + 



32 641/8 
1 1 



V + 



1 



1 1 
81P 64 f/^ ~ 384 ~ 10240 



160 C/5 448 C/7 



(2.42) 



It is noted that the quality factor Q depends on only U. The inner and outer envelope curves are shown by the dotted 
curves in Fig. 3. 

Finally, to use the efficient tiling method developed in [20], we look for an ellipse inside the inner envelope curve 
that can cover the region as much as possible. A simplest choice is to construct the ellipse tangent to the envelope 
curve at dU = with the same curvature. This gives 



2Q4(^4Q2 + i + i) 



v/4Q2 + l(l6Q6 + 8Q4 + Q2 + 
Q^16Q^ + 8Q^ + Q^ + 1) 



(4 02 + 1)2(2 g2 + 1 _ vWTT) 
guu dU^ + gvv dV'^ . 



dU^ 
dV^ 



(2.43) 



This inner ellipse is shown by the thick curve in Fig. 3. It is noted that the ratio between the area covered by the inner 
ellipse and that by the inner envelope curve is 95.0% for Q = 2 and 99.4% for Q = 20, and this ratio is independent 
of the maximum distance dsmax and the central frequency fc- 
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FIG. 3: The dotted curves show the inner and outer envelope curves for the maximum distance ds^ax = 0.02. The thick curve 
is the eUipse given by Eq. (2.43). for all initial phases. Here we set Q = 2. 



E. Conformally flat metric 



Once we have an ellipse that ensures the minimum match we require, it is straightforward to apply our tiling method 
developed in [20]. To do so, wc perform a coordinate transformation that transforms the 2-dimensional metric (2.43) 
to an explicitly conformally flat form. Namely, by the transformation, 



dX = -dV, 
dY = 



9uu 



dU, 



9vv 



we obtain 



dst 



= n{Y) {dx^ + dv^ 



where the conformal factor is given by 

Cl{Y)=gvv{Q{Y)). 
Up to 0(1/ Q^) inclusive, the coordinate transformation is explicitly given as 



X = 
Y = 



1 - 



1 



32 Q« 

1 1 
2Q 



+ 



64 Q8 

3 



+ 



1 



1 



.8Q2 
11 



+ 



64 Q4 

31 



+ 



1 



384 Q6 2048 

69 



16 g2 64 Q3 1024 04 4096 
3357 3891 248531 



32768 Q6 



917504^7 ' 4194304 150994944 Q^' 
To the same accuracy, the inverse transformation becomes 



F = 



i + 2y'' + 3y^ + 



X 



1y^ 

2 



iy3_ly4 
4 32 



ly5_25^6_17^7_ 9447 
4 96 64 7168 



(2.44) 

(2.45) 
(2.46) 



(2.47) 
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O = 1 1 7 9 ^.2 31 97 10649 ^^, 

^ 2y 8 32 128 512 2048 57344 

7457 373751 



■ . (2.48) 



32768 1179648 
The conformal factor ^(y) is given by 

1 3 11 1 ^ 3 ^ . 1 . 98297 

0(F) = ^ H \ \ Y Y^ H Y-^ ^ Y^ 

^ ' 8^2 ley 128 128 2048 4096 229376 

303 ^ 10661897^,^ (2.49) 

224 4718592 ^ ' 

When (5 = 2, the errors induced by the above expansion are found to be ~ 0.1 %. This is accurate enough for our 
purpose as long as we allow the SNR loss, ds^ax> of ^ percent. 

The fact that the conformal factor depends only on the coordinate Y is important. It allows us to apply the same 
tiling method developed in [20], which we recapitulate in Appendix A. In Fig. 8, we show the tiling of the template 
space in the (X, y) coordinates. The tiling in the original coordinates (/c, Q) is shown in Fig. 9. 



III. TEST WITH TAMA NOISE SPECTRUM 



The derivation of the template space metric (2.46) is based on the assumption that the noise is white. This 
assumption may be good because a ringdown wave will be rather narrow banded except for the case Q ~ 2. In order 

to confirm this, we examine the effectiveness of the tiling based on Eq. (2.46) in the case of a colored noise. 

As a model of detector's noise, we use a fitting curve of the one sided noise power spectrum of TAMA300, which is 
given by 

where the frequency / is in units of Hz. This formula of the noise spectrum is obtained during Data Taking 8 in 
2003 [24] . The normalization of the overall amplitude is arbitrary. 

We prepare a template bank using Eq. (2.46). The minimum match is required to be 0.98. We also generate signals 
whose amplitudes arc normalized to unity. Then, we perform the matched filtering and evaluate the maximum of the 
match for each signal. If the match is always greater than 0.98, the use of the template space metric (2.46) is justified 
even for a colored noise, at least for the TAMA DT8 noise spectrum. 

Here we generated 2500 signals randomly with uniform probability in the range of the parameters 1.0 x lO^Hz 
^ /c < 2.5 X lO'^Hz, 2.0 < Q < 33.3 and < 0o < tt. In Fig. 4, we plot the number of signals for each bin of the 
match. We see that the most of the signals are detected with the SNR loss less than 2%. Thus, the template spacing, 
constructed analytically under the assumption of the white noise, turns out to be valid even in the case of the colored 
TAMA noise spectrum. 



IV. DISCUSSION 



The detection of ringdown waves is a direct confirmation of the existence of a black hole. The ringdown waves have 
damped sinusoidal waveforms that reflect the mass and spin of a black hole. In our previous paper [20], we proposed 
an efficient method for tiling the templates for matched filtering in the 2-dimensional {/c, Q} space. However, it 
relied on the template space metric for signals with the initial phase (f>o — 0. In this paper, we took account of the 
initial phase dependence and developed a new, improved method to search for the ringdown waves. 

Since the template metric depends on the initial phase, we first determined the inner envelope curve of all the 
contours of a fixed maximum distance rfs^^x between signals and templates on the (/c, (5)-plane for all possible values 
of the initial phase. We then constructed an ellipse that can cover the region inside the inner envelope as much as 
possible. Finally, we applied our tiling method proposed in [20] to obtain an improved, reliable tiling of the template 
space. 

Another change from our previous work is the difference in the definition of the match. In our previous work, we 
defined the match as the square of the SNR, p^, while here we used the SNR as the match. So, for the same fixed 
maximum distance, the number of templates needed to cover the template space turns out to be smaller than what 
we obtained previously, if we ignore the initial phase dependence. As a result, the actual number of templates to 
cover the template space did not increase much from what was suggested in [20]. 
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0.995 



FIG. 4: The number of signals in terms of the value of the match between the signal and templates. The bank of templates 
are determined assuming the minimum match 0.98 using the metric (2.46) and the tiling method in [20]. The mean value of 
the match is 0.993. 



We also examined the validity of our tiling method in the case of a colored noise spectrum by a Monte Carlo 
simulation. As a model of realistic noise power spectrum, we used a fitting curve of the noise power spectrum of 
TAMA300 during DT8 in 2003. For the pre-assigned maximum allowable SNR loss of 2%, we found that only a few 
out of 2500 signals had the SNR loss larger than 2%. This means that our template spacing is effective even in the 
case of colored noise. 

Since the real data includes the non-stationary, non-Gaussian noise, it is necessary to check the effectiveness of this 
method in a realistic situation. Using the real data of TAMA300, this new template spacing is now being tested by 
Tsunesada et al. [25]. They find a lot of fake events due to the non-stationary, non-Gaussian noise. Particularly, 
the detector has many noise sources that can produce fake ringdown wave signals. Apparently, we need to develop a 
method to remove these fake events without losing real gravitational ringdown wave signals. Perhaps, the best way 
is to perform a coincidence analysis if we have plural detectors. We plan to study methods of coincidence or coherent 
analyses by using several detectors in the future. 
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APPENDIX A: EFFICIENT TILING METHOD 



In Subsec. HE, we have derived the simple, conformally flat metric (2.45) for the template space. Here, using this 
metric, we formulate a tiling algorithm which is not only efficient but also quite simple. 
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1. Basis 



To develop such a method, we note the following. Because of the conformal flatness, the contour of the fixed 

maximum distance cL'P' = (is^j^x c;cntcred at a point on the (X, y)-planc is a circle for sufficiently small ds^-^-^^^. 
Furthermore, along a line of Y =constant, ^{Y) is constant. Thus, choosing first an appropriate Y =constant line, 
say F = gi , we may place circles of the same radius with their centers located along the line F = g'l to cover a region 
surrounding that line. Then, if we find an algorithm to place circles along the Y = qi line and another algorithm 
to choose the next Y =constant line, say Y = q2, to he covered in an appropriate way, we can repeat this tiling 
procedure to cover the whole template space. 

Let us assume that the template space to be tiled is a rectangle given by i^min < F < Fmax and Qmin < Q < Qmax- 
In the {X, Y) coordinates, this rectangle is mapped to the region bounded by the two Y =constant lines corresponding 
to Q = Qmin and Q = Qmax, which we denote hy Y = Yq and Ym, respectively, and the two lines X = Zi-ain{Y) and 
X = Zmax(F) corresponding to = i^min and F = i^max, respectively. Note that Iq > since large Y corresponds 
to small Q. 

First, wc construct a method to determine the spacing of the circles along each Y =constant line. Let us consider 
the line Y = q and place two circles with radius r centered at (p, q) and {p + Ap, q), 

{X-pf + {Y-qf=r\ {X-p-Apf + {Y-qf = r\ (Al) 

We assume Ap < 2r so that the two circles intersect at the two points, (p + Ap/2 , q ± d{r; p)), where d{r; p) is the 
distance to each intersecting point from the line Y = q (see Fig. 5), given by 



d(r; Ap) = yr2-^. (A2) 

Our purpose is to tile the template space by the smallest possible number of filters. In order to do so, we choose the 
parameter p in such a way that the area defined by 5 = Aprf(r; Ap) is maximized, i.e., 

Ap = %/2r, (A3) 
The radius r is determined by the value of ds^^^ and q as 

In this way, we tile the region that covers the line Y = q. 




FIG. 5: Definition of d and r. 



To choose the first line to be covered, we start from the point on the {X, y)-plane corresponding to (F, Q) = 
(-Fmax, Qmin), that Is, {X,Y) = {Xq,Yq) whcrc Xq = ^;max(^)- Then we choose the first line Y = qi and the radius 
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ri so that the point {Xo,Yo) is just on the edge of the first circle and an intersecting point of the first and second 
circles lies on the line Y = Yq as Fig. 6. This is achieved if the center of the first circle is at 



ipu qi) = (Xo - ^l/^/2 , Yq - n/V2^ , 
with the radius determined by solving Eq. (A5), which reads in the present case, 

dsl,,^ = rln{Yo-n/V2). 



(A6) 



(A7) 



0.25- 




FIG. 6: Choosing the first line on the {X, y)-plane. 



Then the center of the n-th circle is at 

iPn, qi) = (Xq - (2n- l)ri/\/2, Fq -ri/\^) , (A8) 
and the number of circles needed to cover the line Y = qi'is given by 

(A9) 

where [x] denotes the maximum integer smaller than a; + 1 and Li is the coordinate length of X to be covered, i.e., 

Ll = Zmax(yo) - -2min(io) 
— -^max -^min 

=: L. (AlO) 

Once the covering of the first line is done, the second Y =constant line is chosen as follows. Let Yi = Y^^ — \/2ri 
and let {Xi,Yi) be the intersecting point of the lines Y = Yi and X = Zmax(F), i-e., {Xi,Yi) = {zniax{Yi),Yi). We 
choose this point as the starting point for the covering of the second line as Fig. 7. That is, the center of the first 
circle {p2, q^) on the second line Y = q2\s 

{P2, q2) = (Xi - r2/\/2 , Fi - r2/V2^ , (All) 

with the radius r2 determined again by Eq. (A5). Then the second line is covered by the same procedure we took for 
the first line. We repeat this procedure until we tile the whole template space we need to cover. 




15 




FIG. 7: Covering of the second line. 



With this tiling procedure, the total mimber of templates is given as follows. Generalizing Eq. (A9), the number of 
templates for the i-th Y =constant line {Y = qi) is 



L 



where is determined by 



ds 



2 2 

— r 



(A12) 



(A13) 



The number of Y =constant lines necessary to cover the template space is determined by the minimum integer v that 
satisfies the inequality, 



(A14) 



i=l 



And the total number of templates is 



(A15) 



i=l 



2. Application 

Let us apply the method developed in the previous subsection to the case of the parameter space (/c, Q) which has 
the range, 

1. X lO^^Hz < fc < 2.5 X lO^Hz, 

2.0 < < 20.0. (A16) 

We set dsf^^^^ = 0.02. This choice is made in order to make the SNR loss to be ^ 3 % in the presence of colored noise 
as discussed in Sec. III. Using Eq.(2.29), we set F ~ ln(/c/100Hz). The above range corresponds to 

3.12 X lO"'' < X < 3.25, 

2.48 X IQ-^ < Y < 2.29 x 10"^ . (A17) 
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if.- Q) 


(,V. V) 


1.0 xlO^Hz,2.0 


3.030840555 xlO'^, 2.293688750 xlO'^ 


1.0 xlO^Hz,20.0 


3.124023844 x 10"", 2.483796010 xlO"^ 


1.0 xlO^Hz,2.0 


2.331909728,2.293688750 xlO"^ 


1.0 xlO^Hz,20.0 


2.302897494, 2.483796010 x 10"^ 


2.5 xlO^Hz,2.0 


3.247808978,2.293688750 xlO'^ 


2.5 xlO^Hz,20.0 


3.219188225,2.483796010 xlO"^ 



TABLE I: The relation between the two coordinates (/c, Q) and {X, Y). 

in the {X, F)-space. (See Table I.) 

Following the procedure described in the previous subsection, we choose the starting point {Xq, Yq) which corre- 
sponds to (/c, Q) = (2.5 X 10'^Hz,2.0). The radius ri is determined by Eq. (A7), or 

om = rln{Yo-n/V2). (Ais) 

Then all the parameters for our tiling method are determined. In Table II, we summarize the tiling parameters. 



i 


(Pi, Qi) 


{Xi, Yi) 


n 









3.247808978, 2.293688750 xlO"^ 






1 


3.202211769, 1.837716656 xlO"^ 


3.229142859, 1.381744561 xlO"^ 


6.448419198 xlO"^ 


36 


2 


3.200527260, 1.095588566 xUr^ 


3.222295278, 8.094325710 xlO"^ 


4.046856890 xlO"^ 


57 


3 


3.205093210, 6.374118861 xlO"^ 


3.219985276, 4.653912012 xlO"^ 


2.432739855 xlO"^ 


94 


4 


3.209939255, 3.649309938 xlO"^ 


3.219230312, 2.644707864 xlO"^ 


1.420721877 xlO"^ 


161 


5 


3.213469023, 2.068578968 xlO'^ 




8.147692989 xlO"^ 


280 



TABLE II: The parameters for the tiling of the template space are summarized. {Xi, Yi) denotes the starting point of the i-th 
template spacing along the line y = Qi. (pi, Qi) denotes the center of the first circle along the i-th line y = qi, ri is the radius 
of the circle, and A'f is the number of circles necessary to cover the i-th line. 

It is noted that the number v of the Y =constant lines is very small, 

!/ = 5. (A19) 

The total number of templates is calculated to be A/^ = 628 {M = 1121 for Qmax = 59.4). In Table III, we summarize 
the total number of templates A/" with the range of the quality factor between Qmin = 2.0 corresponding to a = and 
various maximum Qmax- Here i means the number of Y ^constant lines necessary to cover the template space, and 
the maximum quality factor is determined by this number. 

In Fig. 8, we show the tiling of the template space in the {X, Y) coordinates. The tiling of the template space in 
the original coordinates (/c, Q) is shown in Fig. 9. 

APPENDIX B: PARAMETER ESTIMATION ERRORS 

Here we summarize the method to calculate the parameter estimation errors for gravitational ringdown waves. As 
a signal, we use the ringdown waveforms in Eq. (2.13) with the amplitude A, 



H{A, /e, Q, to, (/-o; /) = A/i((/)o; /c, Q, h; /) , 



(Bl) 



17 



7 




.V 


1 


3.461857533 (0.7045454540) 


36 


2 


6.033964925 (0.9140427132) 


93 


3 


10.60831054 (0.9754675053) 


187 


4 


18.77484416 (0.9931010039) 


348 


5 


33.37367768 (0.9980786204) 


628 


6 


59.48251985 (0.9994680533) 


1121 


7 


106.1826585 (0.9998532378) 


1995 



TABLE III: The maximum range of the quaUty factor and the number of templates. 



0.25 



0.15 




0.05 



FIG. 8: A part of the tiling of the template space in the {X, Y) coordinates. Templates are taken at the centers of the circles. 



and the following calculation is done in the frequency domain. The parameter estimation errors are calculated by 
using Fisher information matrix 



Tjui/ = iH^ij,{A, fc, Q, to, (l)o),H^t^{A, fc, Q, to, ^o)) , 



(B2) 



where /U and v are {InA, to, (po, fc, Q}- In practice, we consider to foto and fc fc/fo with some frequency fo in 
the computational calculation. The root-mean-square errors of parameters are obtained by 



where S = r~^, and the parameter correlation coefhcients c'^'^ between two parameters are 



In order to derive the Fisher information matrix, first we calculate H^ij^ formally as 

1 



(B3) 



(B4) 



-ff,inA(^, fc, Q, to, 4>o; t) = A 



^N{<Po; fc, Q, to) 



Hfc, Q, to, 0o; /) , 
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FIG. 9: A part of the tiling of the template space in the original coordinates (/c, Q) with fc measured in units of lOOHz. Note 
that the contour of ds^ax = 0.02 for each template is warped. 



i?,to(A fc, Q, to, 4>o; t) = A 
H,a{A, fc, Q, to, 0o; t) = A 



1 



1 



iV(0o; fc. Q, to)3/2 
1 

VA^(</.o; fc, Q, to) 



h,toifc, Q, to, 4>o; f) , 

(/l,a(/c, Q, to, 4>o),h{fc, Q, to, (j>o))h{fc, Q, to, (1)0; f) 



h,oc{fc, Q, to, (po; f) 



(B5) 



where a, /3, ... denote {^o, fc, Q}, and we have noted that the normalization constant N does not depend on the 
initial time. And then, the Fisher information matrix is given by 



Tin An — A An ; 



^taW 
ra/3 



1 



7V(0o; fc, Q, to) 



A' 



+ 



N{<l>o; fc, Q, io)2 
1 



{h,to{fc, Q, to, 4>o),h,w{fc, Q, to, 4>o)) ■ 



{h,a{fc, Q, to, 4'o),h{fc, Q, to, <t>o)){hs{fc, Q, to, (t)o),h{fc, Q, to, ^o)) 



N{<l)o; fc, Q, to) 



{h,a,h^l3) 



(B6) 



where W = {to, 0o, fc, Q}- 

Next, we may prepare the partial differentiation h^w of the ringdown waveforms to calculate the above Fisher 
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information matrix. 



Kfc, Q, to, 0o; /) 

h.taifc, Q, to, 0o; /) 
h,4>o{fc, Q, to, (t)o; /) 

hjAhi Q, to, <Po; /) 



h,Q{fc, Q, to, (/'o; /) 



,2 i TT / to 



(cos 00 fc — '^i COS (f)o f Q + 2 fc Q sin (f)o)Qe — ■ 
7T{2fcQ-tfc~2fQ){2f,Q + if, + 2fQ) 



= 2iTT 



fh{.fc, Q, to, 0o; /) , 
2/cQ COS0O - (/c - 2i/<3) sin(/)o ~ ^ . , 
2fcQsm.(j)o + {fc-2ifQ)cos(j)o 



Hfc, Q, to, <Ao; /) 



/ 



(16i//cQ^-4(/, 
-2g(4(/2 + /2)Q2 + /2) sin^o 

(2 (/c + /) Q + i/c) (2 (/c -DQ-ifc) 
((/c-2i/Q) cos^o + 2/c(5 sin(?io) 



X 



/c((2 (/c + /) Q + ^ /c) (2 {fc-f)Q-i fc) cos 0o 



-4i(/,-2i/Q)/,Qsin<^o) 



Hfc, Q, to, (l>o; /) 



(2 (/c + /) + i/c) (2 (/c -f)Q-ifc) 



x{{fc-2ifQ) cos(j)o + 2fcQ sm(j)o)Q 



(B7) 



In practice, by using the inner product (2.7) between two functions, we can obtain the root-mean-square errors and 
correlation coefficients of parameters. 
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